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10 ' Abstract 

00 ■ 

; I We analyse the charge-exchange (CX) measurements 7r~p — )• n^n below pion labo- 

^L) • ratory kinetic energy of 100 MeV. After the removal of five degrees of freedom from 

^vq . the initial database, we combined it with the truncated vr+p database of Ref. [1] and 

fitted to the data using the ETH model [6]. The set of model-parameter values, as 
well as the predictions (which are derived on their basis) for the phase shifts and for 
the low-energy vrA^ constants, are significantly different from the results we obtained 
in the analysis of the truncated n^p elastic-scattering databases. Concerning the 



. - -' phase shifts, the main difference occurs in 6Q_^_ ; smaller differences have been found 

in two p-wave phase shifts. We discuss the implications of these findings in terms of 
the violation of the isospin invariance in the hadronic part of the vrA^ interaction. 
The effect observed is at the level of 7 — 8% in the CX scattering amplitude below 
about 70 MeV. The results and conclusions of this study agree very well with those 
obtained in the mid 1990s, when the isospin invariance was first tested by using ttN 
experimental data, and disagree with the predictions obtained within the framework 
of the heavy-baryon Chiral-Perturbation Theory. 
PACS: 13.75.Gx; 25.80.Dj; 25.80.Gn 
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1 Introduction 



This is the last of three papers addressing issues of the pion-nucleon (ttN) 
interaction at low energies (pion laboratory kinetic energy T < 100 MeV). In 
the first paper [1], we reported on a new phase-shift analysis (PSA) of the 
Ti^p elastic-scattering measurements. In the second paper [2], we examined 
the self-consistency of the large Tr^p elastic-scattering set of Ref. [3], which 
has not been included in our PSAs. In the present study, we will analyse the 
experimental data for the charge-exchange (CX) reaction n^p — ?> n^n and 
investigate whether earlier claims on the isospin breaking [4,5] need, in view 
of the impressive increase of the CX database and of the development of our 
analysis methods during the last fifteen years, to be revised. 

We will mark the physical quantities extracted in the analysis of the ii^p 
elastic-scattering data [1] with the label 'ZUAS12'; this applies both to the 
solution obtained for the parameters of the ETH model [6] and for the hadronic 
phase shifts, as well as to the predictions derived on the basis of Tables 3 (for 
Vmin ~ 1-24 ■ 10~^) and 4 of Ref. [1]. The corresponding results, obtained in 
the present work from the common analysis of the vr+p and CX databases, 
will be marked with the label 'ZUAS12a'. 

Similarly to Ref. [1], we will assume that the physical quantities appearing in 
the present study (i.e., the fit parameters of Sections 4.1 and 4.2, the scatter- 
ing lengths and volumes of Section 4.2.1, the phase shifts of Section 4.2.2, etc.) 
are not purely-hadronic quantities since they still contain residual electromag- 
netic (em) effects. As, however, the repetitive use of the term 'em-modified' is 
clumsy, we will omit it, unless we consider its use necessary as, for instance, 
in the captions of our tables and figures, as well as in Section 7. 



2 Method 



The formalism which we use here has been described in detail in Ref. [7]. 
The determination of the observables from the hadronic phase shifts may be 
found in Section 2 of that work. For n^p scattering, one obtains the partial- 
wave amplitudes from Eq. (1) and determines the no-spin-fiip and spin-flip 
amplitudes via Eqs. (2) and (3). The observables are obtained from these 
amplitudes via Eqs. (13) and (14). For the CX reaction, the observables are 
determined using Eqs. (21-24) of Ref. [7]. 

All the details on the analysis method (i.e., on the minimisation function, on 
the scale factors, etc.) may be found in Section 2.2 of Ref. [1]. The contribution 
xi of the j*^ data set to the overall x^ is given therein by Eq. (1). The scale 



factors Zj, which minimise each x^, are evaluated using Eq. (2); the minimal 
Xj value for each data set (denoted by {x'j)min) is given in Eq. (3) and the 
scaling contribution (of the j*'* data set) to {x'j)min in Eq- (4). Finally, the scale 
factors for free floating Zj (which we will use in Section 5, when investigating 
the absolute normalisation of the CX data using the ZUAS12 prediction as 
reference) are obtained from Eq. (5); their total uncertainty Azj has been 
defined at the end of Section 2.2 of Ref. [1]. 

One statistical test will be performed for each data set, the one involving 
its contribution {x'j)min to the overall x^- The corresponding p-value will be 
evaluated on the basis of {x'j)min and of the number of degrees of freedom of 
the data set (hereafter, the acronym DOF will stand for 'degree(s) of freedom', 
whereas NDF will denote the 'number of DOF'); for a data set with Nj data 
points (none of which is an outlier), NDF is equal to Nj. The p-value for each 
data set will be compared to the confidence level Pmin for the acceptance of 
the null hypothesis (implying no statistically-significant effects). The value 
of Pmm is fixed to the equivalent of a 2.5a effect in the normal distribution, 
corresponding to about 1.24 ■ 10~^. 

The repetitive use of the full description of the databases is largely facilitated 
if we adhere to the following notation: DB+ for the 7i~^p database; DB_ for 
the 7r~p elastic-scattering database; DBq for the CX database; DB_|_/_ for the 
combined vr^p elastic-scattering databases; DB_|_/o for the combined 7r~^p and 
CX databases. 



3 The CX data 



The available CX measurements are listed (in the chronological order in which 
they had been reported) in Table 1. During the last fifteen years, the database 
has been increased sevenfold, i.e., from a mere 47 data points (which were 
available for the analyses [4,5]) to the present status of a total of 333 data 
points. Of the added 286 data points, 270 relate to the differential cross section 
(DCS), 9 to the total cross section (TCS), 6 to the analysing power (AP), and 
1 data point to one threshold constant (to the isovector scattering length bi). 

A milestone in the low-energy CX experimentation was the FITZGERALD86 
[11] experiment which took place at the Clinton P. Anderson Meson Physics 
Facility (LAMPF) and used the 7r° spectrometer to obtain (for the first time) 
important DCS data around the s- and p-wave interference minimum (see next 
section). By establishing a rigorous relation between the real parts of the s- 
and j9-wave CX amplitudes at low energies, the FITZGERALD86 experiment 
became the backbone of the analyses [4,5] and was essential in terms of their 
conclusions on the violation of the isospin invariance in the hadronic part of 



the TrA*" interaction. The FITZGERALD86 experiment was the first complete 
CX experiment, as it also investigated (and reported) the normalisation uncer- 
tainty. Important in terms of the enhancement of the DBq, shortly before the 
millenium change, was the ISENHOWER99 [16] experiment, also performed 
at LAMPF with the 7r° spectrometer. 

Within the last decade, the Crystal-Ball Collaboration made massive contri- 
butions to the low-energy DBq by perfoming experiments at the Brookhaven 
National Laboratory (BNL): the SADLER04 [18] experiment added 60 data 
points, whereas the successor of that experiment, the MEKTEROVIC09 [21] 
experiment, contributed another 140 data points. In total, the data obtained 
by the Crystal-Ball Collaboration amount now to 60% of the entire DBq avail- 
able below 100 MeV. 

From the remaining experiments, the JIA08 [20] data have been taken around 
the s- and p-wave interference minimum, whereas the SCHROEDEROl [17] 
experimental result on the width of the Is state in pionic hydrogen (corrected 
in Ref. [22] after properly taking into account the contributions of the 'yn 
channel) led to the extraction of the scattering length a'^'^ (= \/2bi). The 
FRLEZ98 [14] experiment investigated the angular distribution of the DCS at 
27.50 MeV. Finally, the BREITSCHOPF06 [19] experiment reported the CX 
TCSs at nine energies below 100 MeV. The remaining experiments account 
for less than 10% of the DBq at low energies. 

The full DBq consists of 54 data sets. The quoted values of the TCS in 
Refs. [10,12,21] have been extracted from the coefficients of the Legendre ex- 
pansion of the angular distribution of the DCS and are thus correlated with 
the main results of these experiments; because of this correlation, one may 
use either set of values, but not both. We will use the Legendre coefficients of 
Refs. [10,12] (their DCS measurements have not been published), and directly 
the DCS data of Ref. [21]. 

In our approach, all data sets must be accompanied by a normalisation uncer- 
tainty. This requirement also applies to one-point data sets, because the scale 
factors must be calculated in all cases (in order to enable the investigation of 
a possible bias in the analysis). As a result, we had to assign realistic normal- 
isation uncertainties to those experiments which did not report this quantity. 
We decided to assign these uncertainties as follows: 

• 6% to BUGG71 [8], as the experiments of the TCS or PTCS were assigned 
this uncertainty in Refs. [1,7]. 

• 3% to BREITSCIIOPF06, because the experimental group had already com- 
bined statistical and systematic effects in quadrature (and reported only the 
total uncertainty). 

• 3.1% to SALOMON84 [10], i.e., the normalisation uncertainty of the (very 



similar, as well as close in time) BAGHERI88 [12] experiment. 
• 8% to DUCLOS73 [9], i.e., double the normalisation uncertainty of the 
ISENHOWER99 large-angle data sets. 

The experimental results of Ref. [20] contain asymmetric statistical uncertain- 
ties for the measured DCS. Unable to treat asymmetric uncertainties in our 
analysis, we will use the average of the two uncertainties for each input data 
point. 



4 The s- and p-wave interference minimum 



The main contributions to the CX scattering amplitude in the low-energy re- 
gion come from the real parts of the s and p waves. These contributions are of 
opposite signs and they cancel each other in the forward direction around 45 
MeV. This destructive-interference phenomenon acts as a magnifying glass, 
probing the smaller contributions in the irN dynamics, like those from the 
imaginary parts, from the d and / waves, and (potentially) from isospin break- 
ing. 

An estimation of the energy of the CX DCS minimum has been obtained in 
the FITZGERALD86 experiment, from an extrapolation of their DCS mea- 
surements to centre-of-mass (CM) scattering angle 6 = 0°. According to that 
estimation, the DCS minimum occurs at pion CM kinetic energy e ~ 33 MeV 
or T ~ 45 MeV; the minimal DCS value was also extracted: 2.4±0.5 yub/sr [11]. 
The energy of the s- and p-wave interference minimum has been obtained on 
the basis of the ZUAS12 solution; the prediction is close to 43.7 MeV, whereas 
the minimal DCS falls well below 1 /ib/sr. Both predictions (i.e., of the energy 
of the minimum and of the minimal DCS value) do not seem to match well 
the estimation of Ref. [11]. 

4-1 Fits to the DBq using the K -matrix expansions 

As when analysing the DB_ in Section 3.2 of Ref. [1], we will fix the I = 3/2 
amplitudes from the final fit to the truncated DB_|_ using the i^-matrix expan- 
sions (see Section 3.1 of that paper). The same i^-matrix expansions for the 
1 = 1/2 amplitudes will be used for the description of the CX measurements, 

-1/2 

of course, with different parameters ao+ j ^Ij Ci, c^is, 613, dn, and en. 

The x^ value of the first fit to the data was 400.3 for 326 DOF in the fit, 
indicating a rather coherent database. On the basis of the results of this fit, 
we had to freely float the FITZGERALD86 data set at 40.26 MeV. The second 



run yielded a x^ value of 375.1; we had to freely float the FITZGERALD86 
data set at 36.11 MeV. The third run yielded a x^ value of 351.5; we also had 
to freely float the FITZGERALD86 data set at 32.48 MeV. Interestingly, the 
removal of these three DOF from the initial DBg led to the decrease of the x^ 
by 66.9 units. 

Two additional runs were sufficient in removing all outliers from the DBq. In 
the first run, the one-point data set of the BREITSCHOPF06 experiment at 
75.10 MeV had to be removed, whereas the absolute normalisation of another 
FITZGARALD86 data set (at 47.93 MeV) was removed in the second. The 
final fit to the data yielded the x^ value of 312.8 for 321 DOF in the fit. 
The optimal values of the parameters ao+ , (iis, and eis came out significantly 
different lJ from the ones obtained in the analysis of the DB_ in Ref. [1]. 

Although the results obtained so far provide a firm basis for questioning the 
absolute normalisation in all seven FITZGERALD86 data sets, we will re- 
tain the absolute normalisation of the three remaining FITZGERALD86 data 
sets as their removal is not called for when strictly following our rejection 
criteria. On the other hand, given the importance of FITZGERALD86 data 
in Refs. [4,5], it is obvious that the re-analysis of the CX measurements in 
terms of the violation of the isospin invariance in the hadronic part of the vrA^ 
interaction at low energies is imperative. 



4-2 Common fit to the truncated DB^/q using the ETH model 



Details on the ETH model, as well as on its seven parameters (Gcr, K^, Gp, Kp, 
g-KNN, dnNA, and Z) may be obtained from Refs. [7,1]. This isospin-invariant 
model was introduced in Ref. [6] and was developed to its final form by the 
mid 1990s. 

Prior to fitting the ETH model to the truncated DB_|_/oj we subjected the 
data to a common fit using the i^-matrix expansions. At that moment, the 
truncated DB+ consisted of 340 data points (detailed in Table 1 of Ref. [1]), 
whereas the truncated DBq comprised 328 data points, i.e., the original 333 
data points minus the five outliers detailed in the previous section. The com- 
mon fit to these data using our ii'-matrix expansions resulted in the x^ value 
of 737.0 for 654 DOF in the fit and no outliers. 

The model fit to the data yielded the minimal x^ value of 932.3 for 661 
DOF. The optimal values for the model parameters from the fits to the trun- 
cated DB_|_/o are given in Table 3; the uncertainties contain the Birge factor 



^ The differences between the two sets of values represent effects between 4.9 and 
7.4(T in the normal distribution. 



X^/NDF, which takes account of the goodness of the fit. The table also con- 
tains the ZUAS12 solution for Pmm ~ 1-24 ■ 10~^. The differences between 
these two sets of model-parameter values are significant. 

It is now time to reflect on the important x^ values obtained so far in our 
program. We first concentrate in the results for Pmm ~ 1-24 ■ 10^^. The sep- 
arate fits to the data using the i^-matrix expansions yielded the x^ values of 
427.2, 371.0 [1], and 312.8 (this work) for the truncated DB+, DB_, and DBq, 
with 333, 321, and 321 DOF in the fit, respectively. The x^ values obtained 
with the K-matrix expansions in the two analyses of the combined truncated 
databases (i.e., DB_(_/_ and DB_|_/o) come out very close to the sum of the 
corresponding results for the separate fits: 792.4 (instead of the sum of 798.1) 
for the truncated DB+/_ [1], 737.0 (instead of the sum of 740.0) for the trun- 
cated DB_(./o (this work). (As the NDF in the truncated DB_ and DBq are (by 
chance) identical, the results are directly comparable.) Therefore, we observe 
that, in the case of the fits using the fi'-matrix expansions, the difference of 
the two x^ values (which is about 55.4 units) reflects, almost entirely, the dif- 
ference of the x^ values in the separate fits to the truncated DB„ and DBq 
using the i^- matrix expansions (which is equal to 58.1 units, the smaller x^ 
value corresponding to the fit to the truncated DBq). 

The increase of the x^ values in the fits to the truncated DB_|_/_ using the 
ETH model (over the result of the fits to the same data using the i^-matrix 
expansions) is due to the imposition of theoretical constraints (e.g., of the 
crossing and isospin symmetry); as earlier mentioned, the fits to the data 
using the K-matrix expansions are devoid of theoretical constraints, other 
than the expected low-energy behaviour of the K-matrix elements. All else 
being equal, one would expect that the difference in the x^ values between 
the model fit to the data and the fit using the i^-matrix expansions would be 
(more or less) the same for the two truncated databases, i.e., for the DB_|_/_ 
and DB+/0, inasmuch as the same parameterisation is used in the modelling of 
the two I = 1/2 amplitudes; however, this is far from true. In the case of the 
common fit to the DB_|_/o, the difference between the two x^ values (i.e., for the 
model fit versus the one using i^-matrix expansions) is so large that, though 
the separate fits to the truncated DBq using the i^-matrix expansions yield x^ 
values significantly smaller than those for the truncated DB_, the x^ value for 
the model fit to the truncated DB+/o exceeds by far the one obtained from the 
truncated DB_|_/_. Evidently, the substitution of the truncated DB_ with DBq 
leads to a noticeable deterioration of the overall description of the experimental 
data in the model fits. This deterioration is indicative of a general difficulty in 
the description of the truncated DB+/o in terms of one set of model-parameter 
values. This fact can be explained if the theoretical basis upon which the data 
analysis rests (presumably, the isospin invariance in the hadronic part of the 
ttN interaction) is somewhat disturbed. Inspection of Table 5 shows that the 
information which is obtained from the data for the two other Pmm levels 



used in Ref. [1] (i.e., those corresponding to a 2 and 3cr effect in the normal 
distribution) is very consistent with the resuh for Pmm ~ 1-24 ■ 10~^. We will 
return to this issue in Section 6.3. 



4-2.1 Threshold constants 

From the model parameters and their uncertainties given in Table 3, as well 
as the correlation matrix given in Table 4, we calculated the isoscalar and 
isovector s-wave scattering lengths and the isoscalar(isovector)-scalar(vector) 
p-wave scattering volumes. The results are: 



1 -1/2 , 2 



ao+ + o ao+ = 0.0060(24) fi-\ 



-I «o+ + I «o+ = -0.08272(56) fi-\ 

I -at' + I -at' + I -at' + ^ at' = 0.2085(28) ^-, (1) 

-1 -at' + I at' - I af + j af = 0.1834(17) ^^^ 

^ af + I df - I af - I ~af = -0.1943(17) /x,"^ 

'l af + I at' + I af - ^ af = -0.0680(11) fif. 



Converting these results to the familiar spin-isospin quantities, we obtain 
af = -0.0767(24) fif, hf = 0.1715(28) fif, 

df = 0.2181(19) /i;^ af = -0.0334(10) /i;^ (2) 

a^(2 ^ _o.0443(12) ixf, af = -0.0916(23) fxf. 



Significant differences are found when comparing the values of ag^ , at , and 

— 1/2 

ai_ with the corresponding results of Ref. [1]. 
From the results in Eqs. (2), we obtain 

d^' = ^ -af + i -af = 0.0888(26) ^Jit 



and 



a^° = V2 (-i -al'^ + i ~a?J^^ = -0.11698(79) fi'' 



Unlike the a'^'^ value extracted from the truncated DB+/_ [1], the value of the 
present work is compatible with the result of the measurement of the strong 
shift of the Is state in pionic hydrogen [17]. Additionally, the value of a'^^ is 
marginally consistent (the difference between the two values is at the level of 
l.Scr) with the result of the same experiment for the width of the Is state in 
pionic hydrogen. (In this comparison, we have corrected the raw experimental 
results of Ref. [17] with the em corrections of Ref. [22].) 



4-2.2 Hadronic phase shifts 

The s- and p-wave hadronic phase shifts, from the fit to the truncated DB_|_/o 
using the ETH model, are given in Table 6. These hadronic phase shifts are 
also shown in Figs. 1-6, together with the ZUAS12 results, as well as the cur- 
rent SAID solution (WI08) [23] and their five single-energy values (whenever 
available). A very noticeable difference between the phase-shift results of the 
present work and those obtained in Ref. [1] may be seen in the case of the 5o+ 
phase shift. Smaller differences may be seen in two p-wave phase shifts, in 5i^ 
and 5i_ . 



4-2.3 Scale factors and normalised residuals 

We will next provide evidence that the outcome of the common fit of the 
present work is reliable. Similarly to the tests performed in Section 3.4.4 of 
Ref. [1], we will first investigate whether any bias may be seen in the distri- 
bution of the scale factors Zj, extracted in the final step of the optimisation 
scheme. Subsequently, we will address the issue of the distribution of the nor- 
malised residuals of the fit. 

For both the TT~^p (Fig. 7) and CX (Fig. 8) data sets, the values of Zj which 
lie above and below 1 roughly balance each other and there is no discernible 
energy dependence; a small effect may be present at very low energies (i.e., 
in the scale factors of five of the ISENHOWER99 experiments), yet it does 
not seem to be significant. All taken into account, there is no subrange of the 
entire energy interval in which the data is better or worse fitted than for the 
rest of the range. 

The distribution of the normalised residuals is shown in Fig. 9, along with 
the optimal Gaussian function. The x^ value of this fit was 23.6 for 22 DOF, 



whereas the offset f (for the definition, see Section 3.4.4 of Ref. [1]) comes 
out equal to (—9.3 ±4.0) ■ 10~^. A shift in the distribution, indicating a shght 
underestimation of the 'theoretical' values y'j* in the fit, is visible in Fig. 9. 
However, as the result for f indicates, one may obtain an acceptable description 
of the data by fixing f to 0. For the sake of completeness, we also give the 
optimal value and the uncertainty of parameter B of the Gaussian fit to the 
data: B = 0.501 ± 0.035; the expectation value for B is 0.5. 



Investigation of the absolute normalisation of the CX data on 
the basis of the ZUAS12 solution 



We will next investigate the absolute normalisation of the CX data sets using 
the ZUAS12 solution as reference. To this end, we must determine the amount 
at which the predictions obtained on the basis of the reference solution for each 
CX data set (i.e., the y*^ values appearing in Eq. (1) of Ref. [1]) must be fioated 
in order to optimally fit the data set in question. Therefore, relevant in this 
part of the analysis are the scale factors for free floating Zj, given in Eq. (5) 
of Ref. [1]. 

The extracted values of the scale factors Zj and their total uncertainties may 
be found in Table 7 and, plotted separately for the DCS, TCS, AP, and LEC 
data sets, in Fig. 10. The four FITZGERALD86 data sets which had been 
freely floated in Section 4.1 are not shown; their Zj factors came out equal to 
2.03(14), 2.26(15), 2.28(14), and 1.52(13) (the order corresponds to increasing 
energy); evidently, something must have gone wrong in the determination of 
the absolute normalisation of these data sets. Even when using only the CX 
data, the scale factors obtained for these data sets (see Table 2) are signifl- 
cantly larger than the expectation value of 1. The Zj factor of the one-point 
data set of the BREITSCHOPF06 experiment which was eliminated in Sec- 
tion 4.1 came out equal to 0.927(50), indicating that this measurement is an 
outlier only when compared with the bulk of the CX data; this data point is 
also not shown in Fig. 10. One additional data point, the scattering length 
bi of Ref. [17], is not shown in this flgure. In principle, one could assign this 
data point to the DCS set (after squaring the result for bi), in which case the 
outcome would have been consistent with the scale factors for free floating of 
the ISENHOWER99 10.60 MeV data; however, we show only genuine DCS 
measurements in Fig. 10. 

The careful inspection of Fig. 10 leaves no doubt that, when using the ZUAS12 
solution as baseline, the CX scale factors for free floating contain a large 
amount of fluctuation. Given that the ZUAS12 solution is smooth, the fluctu- 
ation observed in Fig. 10 reflects the scattering of the absolute normalisation 
of the CX data. For instance, the Zj value for the FRLEZ98 data set is equal 



10 



to 1.431(99). This data set lies in between three data sets with significantly 
smaller Zj values, i.e., between the two ISENHOWER99 20.60 MeV data sets 
and the MEKTEROVIC09 33.89 MeV data set. The two neighbouring data 
sets of DUCLOS73 (22.60 and 32.90 MeV), as well as the JIA08 34.37 MeV 
data set have Zj values compatible with 1. 

Restricting ourselves below 70 MeV, we calculated the weighted average of 
the Zj factors of Table 7 and Fig. 10. Our result is that the low-energy CX 
measurements lie on average (15.6 ± 1.4)% above the ZUAS12 predictions. 
Naively translated into a difference in the CX scattering amplitude, this result 
would be equivalent to an effect around the 7 — 8% level. 

In order to provide some perspective and, perhaps, motivation to research 
groups which are active in the low-energy ttN experimentation, we will now 
give the predictions obtained on the basis of the ZUAS12 and ZUAS12a solu- 
tions in a number of situations. We will investigate the differences in these two 
prediction sets and identify the energy regions which provide fertile ground for 
distinguishing experimentally between the two sets of values. We commence 
with the DCSs and TCSs. (The results on the Legendre-expansion coefficients 
are expected to follow the sensitivity of the DCS.) 

The two solutions around the s- and p-wave interference minimum are shown 
in Fig. 11 for 6 = 0°. We observe that the energy corresponding to the DCS 
minimum comes out different in the two solutions; the ZUAS12a prediction 
exceeds the one obtained on the basis of the ZUAS12 results by about 1.2 
MeV. Additionally, the ZUAS12a solution predicts a deeper DCS minimum. 

The shapes of the angular distributions of the CX DCS, obtained on the 
basis of the ZUAS12 and ZUAS12a solutions, are different below and above 
the s- and p-wave interference minimum. Below the minimum, the ZUAS12a- 
based DCS systematically exceeds the one obtained from ZUAS12, by varying 
amounts; at 20 MeV, the relative difference (i.e., the difference normalised to 
the corresponding ZUAS12 values) is 17.7% for ^ = 0° decreasing to 12.2% at 
6 = 180°; the corresponding numbers for 30 MeV are: 24.0 and 11.4%; finally, 
at 40 MeV, the relative differences are: 59.1 and 10.4%. Evidently, the relative 
difference between the two predictions in the forward direction increases as 
the beam energy approaches the energy of the s- and p-wave interference 
minimum. Large effects are also seen in the TCS, slightly decreasing with 
increasing energy, from about 12.9% at 20 MeV, to 11.9% at 30 MeV, and 
to 10.7% at 40 MeV. A representative plot of the two predictions is given in 
Fig. 12 for 30 MeV. 

The picture is slightly different at the higher end of the energies. The ZUAS12a 
result lies here below the ZUAS12 prediction in the forward direction. The 
two predictions cross each other between 40 and 90°; above the crossing, the 
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ZUAS12a prediction exceeds the ZUAS12 one (by smaller amounts when com- 
pared to the low energies). At 60 MeV, the relative difference between the two 
predictions at ^ = 0° is -11.6%, increasing to 8.4% at ^ = 180°; at 80 MeV, 
the values are: —6.0 and 6.4%, whereas at 100 MeV, they are: —5.8 and 4.6%. 
At the same time, the net effect in the TCS decreases with increasing energy, 
from about 7.9% at 60 MeV, to 5.1% at 80 MeV, and to 2.5% at 100 MeV. 
A representative plot of the angular distribution of the CX DCS above the s- 
and j9-wave interference minimum is shown in Fig. 13 for 80 MeV. 

It is evident from our investigation that the AP shows high sensitivity to the 
effect under investigation around the s- and p-wave interference minimum. 
The few available measurements of the AP in the low-energy CX reaction 
have been taken at 98.10 and 100.00 MeV, where the differences between the 
two prediction sets are very small. 

We will now comment on the recent experiment of Ref. [20], which took data at 
forward angles, at six energies around the s- and p-wave interference minimum. 
The authors make the point that their data 'show no evidence for unexpected 
isospin-breaking effects'. To start with, according to Table 2 and assuming 
the correctness of the normalisation of the bulk of the CX data, the JIA08 
measurements lie between 0.5 and 1.50" (o", in this case, is equivalent to 10%) 
below the optimal solution obtained only from the CX data in Section 4.1. We 
created the ZUAS12 and ZUAS12a predictions corresponding to the values of 
the kinematical variables of these data. As seen in Figs. 14, these two solutions 
lie close to one another in the kinematical region of the measurements. As a 
result, the experiment indeed agrees with the ZUAS12 solution, but it also 
does with the ZUAS12a prediction (which was not available at the time the 
experimental paper appeared). Given the large normalisation uncertainties of 
the low-energy irN experiments, as well as the general closeness of the ZUAS12 
and ZUAS12a solutions, it is rather unlikely that any single experiment will 
be capable of disproving the violation of the isospin invariance in the hadronic 
part of the ttN interaction, especially one with a normalisation uncertainty at 
the 10% level. 

We would like to stress that our predictions for the usual observables (DCS, 
TCS, AP), calculated on the basis of the ZUAS12 and ZUAS12a solutions, 
are available to all researchers in the iiN field upon request; it takes very 
little effort to produce reliable results from our software on a grid of {0,T) 
points. We must, however, draw attention to the correct interpretation of 
our numbers and figures both of the present work, as well as of any future 
communication (to irN research groups) as a response to a specific request. 
An agreement of the results of a CX experiment with the ZUAS12a solution 
may not be interpretted as unmistakable evidence for the violation of the 
isospin invariance in the hadronic part of the vrA^ interaction (though such 
an agreement is in the right direction); it simply shows compatibility of the 
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experiment in question with the bulk of the CX data as analysed here. On the 
other hand, a significant disagreement of a CX experiment with the ZUAS12a 
solution points in the direction of potential problems with the integrity of the 
CX data; if, in this last case, the experimental results agree with the ZUAS12 
prediction, then the interpretation can only be that the experiment does not 
support the possibility of the violation of the isospin invariance in the hadronic 
part of the vrA^ interaction as it emerges in the present study. 



Possible causes of the observed differences between the ZUAS12 
and ZUAS12a solutions 



We now summarise the main results obtained so far. 

• Three of the values of the parameters of the A'-matrix expansions (i.e., 
ao+ ! di3, and Cis), extracted in the analysis of the truncated DB+/01 differ 
significantly from those obtained in the analysis of the truncated DB+/_. 

• Two of the values of the parameters of the ETH model (i.e., Gp and QnNN), 
extracted in the analysis of the truncated DB+/O) differ significantly from 
those obtained in the analysis of the truncated DB+/_. 

• When using the ETH model in the common fits to the data, the substitution 
of the truncated DB_ with DBq leads to obvious deterioration of the results, 
indicating difficulties in the description of these measurements in terms of 
one set of model-parameter values. In this respect, the results of Table 5, for 
the three values of Pmin which have been used in Ref. [1], are very consistent. 

• A significant difference has been seen in the s-wave phase shift 6q[^ ; milder 

~l/2 

differences are observed in two p-wave phase shifts, i.e., in the 6i_ and in 
the 5i('. 

• The reproduction of the absolute normalisation of the CX data sets on 
the basis of the ZUAS12 solution is very poor below about 70 MeV. The 
differences amount to a 7 — 8% effect in the CX scattering amplitude. 

These differences between the ZUAS12 and ZUAS12a solutions and the predic- 
tions obtained on their basis are significant. We will now attempt to identify 
possible causes of the observed discrepancies. 

In view of these results, there are three hypotheses which cannot all be simul- 
taneously vahd. 

• The bulk of the low-energy iiN experimental data is reliable. 

• The missing parts in the em corrections, applied to the experimental data 
in order to extract the hadronic part of the iiN amplitude, are negligible. 

• The isospin invariance in the hadronic part of the ttN interaction is obeyed. 
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We will now elaborate on these three possibilities. 
6.1 Experimental problems 



The first explanation for the observed differences involves a trivial effect, i.e., 
the systematic incorrectness (overestimation or underestimation) of the abso- 
lute normalisation in the low-energy vrA^ data. Furthermore, arguments have 
been presented in Ref. [7], to support the claim that the reported experimen- 
tal uncertainties are underestimated. In this respect, the two elastic-scattering 
reactions seem to be more affected (than the CX). 

Our first point concerns the statistical uncertainties of the data points. When 
visually inspecting the low-energy ttN data, one is frequently unable to com- 
prehend how it is possible for successive measurements (i.e., at neighbouring 
values of the CM scattering angle) to be so different. There are many occasions 
in the databases when one starts questioning the statistical uncertainties of the 
data points, which, in fact, should be the easiest to obtain in an experiment. 

Our second point concerns the systematic effects. It is not understood how 
the absolute normalisation of some experiments, e.g., of the FITZGERALD86 
data sets at the three lowest energies, may be wrong (according to the bulk of 
the truncated DBq) by an average of about 70% (and of a fourth one by 45%), 
at a time when the reported normalisation uncertainty in the experiment was 
7.8%. Such an effect may only be caused because of any of three reasons (or 
their combination): a) the determination of the absolute normalisation in the 
experiment has been erroneous, b) the energy of the incoming beam has not 
been what the experimenters expected, or c) the normalisation uncertainty in 
the experiment has been grossly underestimated. 

Our third point concerns the small values of the normalisation uncertainty 
reported in many low-energy vrA^ experiments; for instance, the reported nor- 
malisation uncertainty in 37 out of the 90 data sets in the initial DB+/_ is 
below 3%. One could possibly try to manipulate the normalisation uncertainty 
reported in such experiments (perhaps, by setting these values to 3%), yet such 
an approach would be arbitrary. Therefore, one must reluctantly accept the 
reported values. As a result, one must rely on an approach which places impor- 
tance on the bulk of the experimental data and apply a reasonable procedure 
for the elimination of the outliers; this is the approach which we set up in 
Refs. [1,2,7]. 

There is one disturbing discrepancy in the entire analysis of the truncated 
DB_|_/_ which cannot be easily put aside, i.e., the disagreement between the 
a'^'^ value obtained as an extrapolation from the experimental data (above 
the TiN threshold) and the one extracted directly at the nN threshold from 



14 



the strong shift of the Is level in pionic hydrogen. Assuming the correctness 
of both the bulk of the truncated DB_ and of the raw measurement of ei^ 
[17], the two extracted values should be compatible, if a consistent set of 
em corrections (i.e., corrections which have been obtained within the same 
framework), which is also complete in the sense of containing the contributions 
from all relevant physical effects, have been applied to the raw data. In case 
of any missing pieces in the em corrections (see next section), the question 
surfaces as to the energy dependence of the missing contributions. At present, 
there are no known effects which could invalidate the result from the pionic- 
hydrogen experiment [17] and, presumably, it is equally difficult to accept that 
the bulk of the truncated DB_ is not reliable, especially as this set comprises 
measurements which 'cluster' together better than the ones for 7r+p elastic 
scattering. 

On the other hand, though there might be problems with the bulk of the low- 
energy TiN experimental data, it is not at all clear how they can create the 
results of Table 5. 



6.2 Missing pieces in the em corrections 



Although it is not clear how these contributions can modify so drastically the 
overall picture and especially the results of Table 5 (any missing em contri- 
butions should be common in the description of the experimental data on the 
basis of the i^-matrix expansions and of the ETH model), this is an impor- 
tant issue which must be properly defined and addressed. In Refs. [7,22], some 
details are given on the effects which the stage-II em corrections may contain. 
These effects are mostly related to the use of the (unknown) hadronic masses 
for the proton, the neutron, and the charged and neutral pion in the determi- 
nation of the em corrections; the physical masses have been used in Refs. [28] , 
which we apply in our program. Lacking, however, a consistent method for 
obtaining even a rough estimation of these contributions makes this discus- 
sion rhetorical and speculative. Until these values become available, there is 
no other option but to assume their smallness. 



6.3 The violation of the isospin invariance in the hadronic part of the tcN 
interaction 



This is the third of the available options and, admittedly, the most exciting one 
in Physics terms. Interestingly, this possibility may easily provide an explana- 
tion for the results of Table 5. There has been a great amount of discussion 
regarding the acceptance of the conclusions of Refs. [4,5]. One, however, is 
tempted to pose the question: 'Why should the isospin invariance in the 7iN 
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system be obeyed in the first place?' After all, the hadronic masses of the u 
and d quarks are different; similarly, the masses of the nucleons differ (beyond 
trivial em effects), and so do those of A's. It appears, therefore, that the proper 
question is not whether the isospin invariance is violated, but at which amount 
it is. Within the framework of the heavy-baryon Chiral-Perturbation Theory, 
the group of Meif5ner has repeatedly treated isospin-breaking effects in the iiN 
system at low energies and invariably concluded that such effects should be 
small, at the percent level at most (e.g., see Ref. [24] and the references cited 
therein) . 

Two mechanisms had been proposed in the past, to account for the violation 
of the isospin invariance in the hadronic part of the irN interaction at the level 
of Feynman diagrams: the first mechanism affects the elastic scattering {p^ — u 
mixing [25]), the second the CX reaction {rj—ir^ mixing [26]). As both the u and 
the f] states are singlets, the coupling of the former with the p" and of the latter 
with the 7r° explicitly violate the isospin invariance. Given that, in the case 
of the elastic scattering, only one graph (i.e., the t-channel p-exchange graph) 
is affected, whereas in the case of the CX reaction all Feynman diagrams are 
affected (see Fig. 15), one would be tempted to jump to the conclusion that the 
isospin-breaking effects are expected to be more important in the latter case; 
however, Ref. [26] concluded with the statement that '. . .the isospin violation 
from f] — 71^ mixing can be safely ignored in vrA^ partial- wave analyses.' 

We noticed that the coupling constant q^nn is significantly affected when sub- 
stituting the truncated DB_ with DBq. Of course, if the isospin invariance is 
violated, there is not one coupling constant Qt^nn] one must distinguish be- 
tween g^^Opp, gjrOnn, and g.„±pn- In this case, the fits to the elastic-scattering data 
essentially determine g-,^±pn-, whereas those involving the CX reaction will also 
contain contributions from g-j^Opp and gnOnn- As a result, the value of the cou- 
pling constant g^NN extracted from the common fits to the truncated DB_|_/o 
represents a weighted average of these three quantities. The differences ob- 
served imply that at least one of the two g-^ONN coupling constants is different 
from g-jripn- The isospin-breaking effects on these coupling constants have been 
studied theoretically and generally found small; for instance, Ref. [27] evalu- 
ated these effects with QCD sum rules and reported that g-jripn should be 
equal to the average of the two g^oNN values, and that the splitting should be 
expected between 1.2 and 3.7%. The difference between the two g-^NN values 
of Table 3 is larger, around the 4.5% level. 



7 Discussion 



The present study concludes our program of re-analysing the pion-nucleon 
{t[N) data below pion laboratory kinetic energy of 100 MeV. The separate 
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analysis of the data for the two elastic-scattering and for the charge-exchange 
(CX) reactions was enabled via suitable expansions of the s- and p-wave K- 
matrix elements at low energies. Common fits to the data were performed 
using either these i^-matrix expansions or the ETH model of Ref. [6], which is 
based on meson t-channel exchanges, as well as on the s- and -u-channel N and 
A contributions. As a principle, the analysis with the i^-matrix expansions 
was used in order to identify the obvious outliers in the databases and to test 
the self-consistency of the input prior to its submission to the fits using the 
ETH model. Given the freedom rendered to the input data in the analysis 
with the i^-matrix expansions, any outliers in the fits can only be attributed 
to experimental problems. After the removal of the outliers, the databases 
were submitted to common fits using the ETH model. The optimal values of 
the model parameters, as well as their correlation matrix, are obtained from 
these fits and can be used as the basis for generating Monte-Carlo predictions 
for the phase shifts, for the low-energy nN constants, as well as for the stan- 
dard observables (i.e., for the differential cross section (DCS), analysing power 
(AP), partial-total cross section (PTCS), and total cross section (TCS)) for 
any of the three reactions, at any value of the relevant kinematical variables 
(i.e., energy and scattering angle for the DCS and AP, energy and angle limit 
for the PTCS, and energy for the TCS). 

Given that the electromagnetic (em) corrections (which are applied to the 
hadronic phase shifts in order to extract the vrA^ partial-wave amplitudes, 
which, in turn, lead to the observables) of Ref. [28] have been obtained by 
using the physical, instead of the (unknown) hadronic, masses for the proton, 
the neutron, and the charged and neutral pion, we have assumed the cautious 
attitude of considering our physical quantities (i.e., the model parameters, the 
low-energy vrA^ constants, the phase shifts, etc.) not purely hadronic, but em- 
modified hadronic. At the present time, one cannot give an estimation of the 
importance of the missing em effects (i.e., of the stage-H em corrections). 

Following the procedure described in the first paragraph of the present section, 
we first analysed the two elastic-scattering reactions and obtained the solution 
(for the model parameters, as well as for the em-modified phase shifts) [1] . In 
this paper, we analysed the n^p and CX databases. By comparing the results 
of these two PSAs, we found large effects both in the model parameters, as well 
as in the em-modified hadronic phase shifts. We have provided evidence that 
the differences are predominantly attributable to the s waves, in particular, to 
the em-modified phase shift (5o+ ; smaller differences have been found in two 
p-wave em-modified phase shifts. 

Assuming the integrity of the bulk of the low-energy ttN databases, as well as 
the smallness of the missing pieces in the em corrections (i.e., of the stage-H 
effects), these discrepancies can only be blamed on the violation of the isospin 
invariance in the hadronic part of the nN interaction at low energies. The 
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effect observed is at the level of 7 — 8% in the CX scattering amplitude below 
about 70 MeV. 

The results of this study agree very well with those obtained in the mid 1990s, 
when the isospin invariance in the ttN system was first tested by using the 
available experimental information. This agreement is remarkable given the 
changes of the databases in the meantime (e.g., the CX database has been 
increased sevenfold), the analysis methods, and the em corrections which are 
applied to the experimental input data. Our results are in severe disagreement 
with predictions (for the isospin-breaking effect in the CX scattering ampli- 
tude), obtained within the framework of the heavy-baryon Chiral- Perturbation 
Theory; according to these estimations, the expected effects should be around 
the percent level. 

In any analysis of experimental data, a question surfaces sooner or later: to 
what extent can the bulk of the data be trusted? The answer to this question 
in the case of the low-energy vrA^ interaction is not simple. For a few years, we 
believed that, by excluding a small number of obvious outliers, we could obtain 
self-consistent databases. However, the analysis of the recently-obtained Denz 
et al. [3] measurements, which comprise a database of their own, led to numer- 
ous problems [2] mainly relating to the shape of the angular distributions of 
their ir^p differential cross sections. From the experimental point of view, it is 
as if we have not really advanced much after three decades of experimentation. 
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Table 1 

The low-energy CX experiments in chronological order. The first column contains a 
label identifying the experiment. Columns 2—5 contain the number of measurements 
reported in each experiment: DCS stands for the differential cross section, LEC for 
the first three coefficients in the Legendre expansion of the DCS, TCS for the total 
cross section, and AP for the analysing power. The adjacent column contains the 
energy or energy range of the experiment; the centre-of-mass scattering angle or 
angular interval of the measurements is listed in the last column. The experiment 
of Ref. [17] extracted the scattering length a'^^ (= v2&i) from a measurement of 
the width of the Is state in pionic hydrogen; thus, it cannot be placed in any of the 
categories under Columns 2 — 5. 



Experiment 


DCS 


LEC 


TCS 


AP 


T (MeV) 


e 


BUGG71 [8] 






1 




90.90 




DUCLOS73 [9] 


3 








22.60 - 42.60 


180° 


SALOMON84 [10] 




6 






27.40, 39.30 




FITZGERALD86 [11] 


21 








32.48 - 63.21 


9.60 - 25.04° 


BAGHERI88 [12] 




12 






45.60-91.70 




STASK093 [13] 








4 


100.00 


75.00 - 130.00° 


FRLEZ98 [14] 


6 








27.50 


4.70 - 50.90° 


GAULARD99 [15] 








6 


98.10 


8.02 - 86.05° 


ISENHOWER99 [16] 


40 








10.60 - 39.40 


9.60 - 168.24° 


SCHROEDEROl [17] 










0.00 




SADLER04 [18] 


60 








63.86 - 94.57 


18.19- 161.81° 


BREITSCHOPF06 [19] 






9 




38.90 - 96.50 




JIA08 [20] 


24 








34.37 - 59.68 


5.81-41.39° 


MEKTEROVIC09 [21] 


140 








33.89 - 86.62 


18.19- 161.81° 
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Table 2 

The data sets comprising the truncated database for the CX reaction, the pion 
laboratory kinetic energy T (in MeV), the number of degrees of freedom {NDF)j 
for each data set, the scale factor Zj which minimises x? (Eq. (1) of Ref. [1]), the 
values of {x'j)min, and the p- value of the fit for each data set. The numbers of this 
table correspond to the final fit to the data using the i^-matrix expansions (see 
Section 4.1). 



Data set 


T 


{NDF)j 


Zj 


{Xj Jmin 


p- value 


Comments 


BUGG71 


90.90 


1 


1.0225 


0.1470 


0.7014 




DUCLOS73 


22.60 


1 


0.9413 


1.2465 


0.2642 




DUCLOS73 


32.90 


1 


0.9717 


0.2700 


0.6033 




DUCLOS73 


42.60 


1 


0.9098 


2.3836 


0.1226 




SALOMON84 


27.40 


3 


0.9720 


2.8685 


0.4124 




SALOMON84 


39.30 


3 


0.9937 


1.0774 


0.7825 




FITZGERALD86 


32.48 


2 


1.5076 


2.3635 


0.3067 


freely floated 


FITZGERALD86 


36.11 


2 


1.7103 


1.1845 


0.5531 


freely floated 


FITZGERALD86 


40.26 


2 


1.8274 


6.1362 


0.0465 


freely floated 


FITZGERALD86 


47.93 


2 


1.4497 


1.5402 


0.4630 


freely floated 


FITZGERALD86 


51.78 


3 


1.1236 


7.3728 


0.0609 




FITZGERALD86 


55.58 


3 


1.0926 


2.5611 


0.4644 




FITZGERALD86 


63.21 


3 


1.0503 


1.2246 


0.7471 




BAGHERI88 


45.60 


3 


1.0056 


0.1314 


0.9878 




BAGHERI88 


62.20 


3 


0.9589 


3.4999 


0.3208 




BAGHERI88 


76.40 


3 


0.9731 


3.2706 


0.3518 




BAGHERI88 


91.70 


3 


1.0151 


2.8032 


0.4230 




STASK093 


100.00 


4 


0.9948 


1.4336 


0.8383 




FRLEZ98 


27.50 


6 


1.0902 


10.2313 


0.1152 




GAULARD99 


98.10 


6 


1.0241 


1.1007 


0.9815 




ISENHOWER99 


10.60 


4 


1.0203 


2.1816 


0.7024 




ISENHOWER99 


10.60 


5 


1.0054 


1.4611 


0.9175 




ISENHOWER99 


10.60 


6 


1.0181 


8.0844 


0.2320 




ISENHOWER99 


20.60 


5 


0.9803 


1.5435 


0.9080 




ISENHOWER99 


20.60 


6 


1.0120 


8.1813 


0.2251 
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Table 1 continued 



Data set 


T 


{NDF)j 


^j 


\Xj )min 


p- value 


Comments 


ISENHOWER99 


39.40 


4 


1.0701 


7.1132 


0.1300 




ISENHOWER99 


39.40 


5 


1.0597 


8.4184 


0.1346 




ISENHOWER99 


39.40 


5 


0.9514 


5.1617 


0.3965 




SCHROEDEROl 


0.00 


1 


0.9747 


2.5184 


0.1125 




SADLER04 


63.86 


20 


0.9548 


16.0803 


0.7116 




SADLER04 


83.49 


20 


0.9881 


11.6506 


0.9276 




SADLER04 


94.57 


20 


1.0296 


7.2573 


0.9958 




BREITSCHOPF06 


38.90 




0.9960 


0.1643 


0.6852 




BREITSCHOPF06 


43.00 




1.0011 


0.0259 


0.8721 




BREITSCHOPF06 


47.10 




0.9981 


0.0572 


0.8110 




BREITSCHOPF06 


55.60 




0.9952 


0.2074 


0.6488 




BREITSCHOPF06 


64.30 




0.9725 


3.7739 


0.0521 




BREITSCHOPF06 


65.90 




0.9779 


2.3441 


0.1258 




BREITSCHOPF06 


76.10 




0.9814 


1.6114 


0.2043 




BREITSCHOPF06 


96.50 




0.9816 


0.6152 


0.4328 




JIA08 


34.37 




0.8434 


4.9306 


0.2945 




JIA08 


39.95 


4 


0.8680 


3.1715 


0.5295 




JIA08 


43.39 


4 


0.8777 


2.5167 


0.6416 




JIA08 


46.99 


4 


0.9798 


5.1175 


0.2754 




JIA08 


54.19 


4 


0.9080 


2.0430 


0.7279 




JIA08 


59.68 


4 


0.9265 


3.2449 


0.5177 




MEKTEROVIC09 


33.89 


20 


1.0239 


17.0075 


0.6525 




MEKTEROVIC09 


39.38 


20 


1.0145 


14.7514 


0.7905 




MEKTEROVIC09 


44.49 


20 


1.0100 


33.1457 


0.0325 




MEKTEROVIC09 


51.16 


20 


1.0357 


15.0473 


0.7737 




MEKTEROVIC09 


57.41 


20 


1.0394 


19.9034 


0.4640 




MEKTEROVIC09 


66.79 


20 


1.0235 


19.4707 


0.4914 




MEKTEROVIC09 


86.62 


20 


1.0019 


31.1877 


0.0528 
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Table 3 

The values of the seven parameters of the ETH model obtained from the common fit 
to the truncated Tr~^p and CX databases (ZUAS12a solution). The ZUAS12 solution 
[1], obtained from the fits to the it p elastic-scattering data for Pmin ~ 1-24 • 10~^, 
is shown for comparison. 





This work (ZUAS12a) 


ZUAS12 


G^iGeV-^) 


29.6 ±1.9 


27.48 ± 0.86 


K. 


0.117 ±0.054 


0.016 ± 0.034 


GpiGeV-^) 


59.36 ± 0.57 


54.67 ±0.61 


Kp 


1.34 ±0.29 


0.66 ± 0.41 


g-n-NN 


13.43 ±0.11 


12.84 ±0.12 


g-KNA 


29.07 ±0.28 


29.77 ±0.26 


z 


-0.393 ± 0.086 


-0.552 ±0.056 



Table 4 

The correlation matrix for the seven parameters of the ETH model for the common 
fit to the truncated n'^p and CX databases. 





G„ 


K, 


G, 


Kp 


g-KNN 


g-n-NA 


Z 


G. 


1.0000 


0.5620 


0.1167 


0.2311 


0.3777 


-0.4932 


-0.2225 


K^ 


0.5620 


1.0000 


0.7092 


0.8642 


0.8280 


-0.9501 


0.6588 


Gp 


0.1167 


0.7092 


1.0000 


0.7219 


0.8386 


-0.7537 


0.7780 


Kp 


0.2311 


0.8642 


0.7219 


1.0000 


0.7628 


-0.8591 


0.8060 


g-TTNN 


0.3777 


0.8280 


0.8386 


0.7628 


1.0000 


-0.9024 


0.6317 


gwNA 


-0.4932 


-0.9501 


-0.7537 


-0.8591 


-0.9024 


1.0000 


-0.6497 


z 


-0.2225 


0.6588 


0.7780 


0.8060 


0.6317 


-0.6497 


1.0000 
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Table 5 

The various x^ values obtained in the analysis of the low-energy vrAT data, along with 
the number of degrees of freedom in each fit, for three values of Pmin (the confidence 
level in the statistical tests); these three Pmin values correspond to a 3, 2.5, and 2a 
effect in the normal distribution, respectively. The definitions of the databases (DB) 
are given at the end of Section 2. Separate fits to the three databases using the ETH 
model are not attempted, given the largeness of the correlations among the model 
parameters in that case. 



Parametric model 


DB+ [1] 


DB_ [1] 


DBo (this work) 


DB+/_ [1] 


DB_|_/o (this work) 


Prmn ~ 2.70 • IQ-^ (Scj) 


iT-matrix 
ETH model 


434.4/334 


397.9/325 


333.4/323 


825.9/659 
905.4/666 


765.1/657 
969.8/664 


Pmin ~ 1.24 • 10-2 (2.5ct) 


i^-matrix 
ETH model 


427.2/333 


371.0/321 


312.8/321 


792.4/654 
872.9/661 


737.0/654 
932.3/661 


Prmn ~ 4.55 • 10-^ {2a) 


iT-matrix 
ETH model 


357.0/310 


332.3/316 


306.4/320 


684.4/626 
755.2/633 


663.3/630 
814.6/637 
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Table 6 

The values of the six s- and p-wave em- modified hadronic phase shifts (in degrees), 
obtained on the basis of the results of Tables 3 (ZUAS12a solution) and 4. 



T{MeV) 


6li' (S31) 


4'+ (sii) 


~5li' (P33) 


6f_^ (P31) 


~6li' (P13) 


Sl^_^ (PU) 


20 


-2.456(49) 


4.565(61) 


1.3239(89) 


-0.2390(71) 


-0.1694(62) 


-0.430(13) 


25 


-2.861(51) 


5.092(65) 


1.877(12) 


-0.329(10) 


-0.2293(87) 


-0.569(18) 


30 


-3.260(52) 


5.562(68) 


2.507(15) 


-0.425(13) 


-0.292(11) 


-0.708(24) 


35 


-3.658(53) 


5.988(71) 


3.212(17) 


-0.528(17) 


-0.357(15) 


-0.843(30) 


40 


-4.055(53) 


6.376(73) 


3.996(20) 


-0.635(21) 


-0.423(18) 


-0.973(37) 


45 


-4.455(54) 


6.733(76) 


4.861(22) 


-0.747(25) 


-0.489(22) 


-1.095(44) 


50 


-4.858(56) 


7.062(81) 


5.810(24) 


-0.863(30) 


-0.556(25) 


-1.208(52) 


55 


-5.264(58) 


7.366(86) 


6.849(25) 


-0.982(35) 


-0.622(29) 


-1.310(60) 


60 


-5.673(62) 


7.647(92) 


7.985(27) 


-1.104(40) 


-0.688(34) 


-1.401(69) 


65 


-6.087(68) 


7.907(99) 


9.222(30) 


-1.229(46) 


-0.753(38) 


-1.479(78) 


70 


-6.504(75) 


8.15(11) 


10.571(33) 


-1.356(52) 


-0.817(43) 


-1.543(87) 


75 


-6.925(84) 


8.37(12) 


12.038(38) 


-1.486(58) 


-0.880(48) 


-1.593(98) 


80 


-7.350(94) 


8.57(13) 


13.635(45) 


-1.617(64) 


-0.942(54) 


-1.63(11) 


85 


-7.78(11) 


8.76(14) 


15.372(55) 


-1.751(71) 


-1.003(59) 


-1.65(12) 


90 


-8.21(12) 


8.93(16) 


17.261(68) 


-1.885(79) 


-1.061(65) 


-1.65(13) 


95 


-8.65(13) 


9.09(17) 


19.315(84) 


-2.022(86) 


-1.118(71) 


-1.64(14) 


100 


-9.08(15) 


9.23(19) 


21.55(10) 


-2.159(94) 


-1.174(77) 


-1.61(16) 
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Table 7 

The scale factors Zj for free floating for the CX data sets, obtained on the basis of the 
ZUAS12 solution [1], listed separately for the differential cross sections (upper part), 
total cross sections (second part), analysing powers (third part), and the results for 
the coefficients of the Legendre expansion of the differential cross section (last part). 
The four FITZGERALD86 data sets which had been floated freely in Section 4.1 
have not been included. Not listed in the table is also the result of Ref. [17] for the 
isovector scattering length bi. The quantity Azj is the total uncertainty (see end of 
Section 2.2 of Ref. [1]). 



Data set 


T (MeV) 


% 


Azj 


DUCLOS73 


22.60 


1.03 


0.14 


DUCLOS73 


32.90 


1.09 


0.13 


DUCLOS73 


42.60 


0.95 


0.12 


FITZGERALD86 


51.78 


1.29 


0.10 


FITZGERALD86 


55.58 


1.252 


0.099 


FITZGERALD86 


63.21 


1.204 


0.093 


FRLEZ98 


27.50 


1.431 


0.099 


ISENHOWER99 


10.60 


1.45 


0.12 


ISENHOWER99 


10.60 


1.332 


0.081 


ISENHOWER99 


10.60 


1.307 


0.057 


ISENHOWER99 


20.60 


1.205 


0.052 


ISENHOWER99 


20.60 


1.229 


0.047 


ISENHOWER99 


39.40 


1.47 


0.11 


ISENHOWER99 


39.40 


1.245 


0.045 


ISENHOWER99 


39.40 


1.089 


0.042 


SADLER04 


63.86 


1.050 


0.068 


SADLER04 


83.49 


1.045 


0.053 


SADLER04 


94.57 


1.063 


0.047 


JIA08 


34.37 


1.04 


0.11 


JIA08 


39.95 


1.00 


0.12 


JIA08 


43.39 


0.94 


0.13 


JIA08 


46.99 


1.09 


0.14 


JIA08 


54.19 


0.94 


0.13 


JIA08 


59.68 


0.99 


0.12 
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Table 6 continued 



Data set 


T (MeV) 


% 


Azj 


MEKTEROVIC09 


33.89 


1.213 


0.040 


MEKTEROVIC09 


39.38 


1.182 


0.032 


MEKTEROVIC09 


44.49 


1.160 


0.031 


MEKTEROVIC09 


51.16 


1.181 


0.033 


MEKTEROVIC09 


57.41 


1.168 


0.031 


MEKTEROVIC09 


66.79 


1.127 


0.032 


MEKTEROVIC09 


86.62 


1.056 


0.030 


BUGG71 


90.90 


1.068 


0.061 


BREITSCHOPF06 


38.90 


1.12 


0.10 


BREITSCHOPF06 


43.00 


1.17 


0.15 


BREITSCHOPF06 


47.10 


1.11 


0.12 


BREITSCHOPF06 


55.60 


1.078 


0.094 


BREITSCHOPF06 


64.30 


0.967 


0.069 


BREITSCHOPF06 


65.90 


0.995 


0.067 


BREITSCHOPF06 


76.10 


0.994 


0.065 


BREITSCHOPF06 


96.50 


0.999 


0.039 


STASK093 


100.00 


0.91 


0.11 


GAULARD99 


98.10 


0.961 


0.058 


SALOMON84 


27.40 


1.099 


0.056 


SALOMON84 


39.30 


1.132 


0.059 


BAGHERI88 


45.60 


1.149 


0.036 


BAGHERI88 


62.20 


1.042 


0.039 


BAGHERI88 


76.40 


1.038 


0.036 


BAGHERI88 


91.70 


1.063 


0.039 
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Fig. 1. The em-modified hadronic phase shift 6q_j_ (S31) from the present work, 
along with la uncertainties (yellow band). Shown is also the ZUAS12 solution [1], 
obtained on the basis of the elastic-scattering data below 100 MeV, along with the 
corresponding la uncertainties (blue band). The current SAID solution (WI08) [23] 
is represented by the dashed curve; their five single-energy values (at T = 20, 30, 
47, 66, and 90 MeV) are also shown. 
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~l/2 

Fig. 2. The em-modified hadronic phase shift (5q1j_ (Sll) from the present work, 
along with la uncertainties (yellow band). Shown is also the ZUAS12 solution [1], 
obtained on the basis of the elastic-scattering data below 100 MeV, along with the 
corresponding la uncertainties (blue band). The current SAID solution (WI08) [23] 
is represented by the dashed curve; their five single-energy values (at T = 20, 30, 
47, 66, and 90 MeV) are also shown. 
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Fig. 3. The em-modified hadronic phase shift 5i_^ (P33) from the present work, 
along with la uncertainties (yellow band). Shown is also the ZUAS12 solution [1], 
obtained on the basis of the elastic-scattering data below 100 MeV, along with the 
corresponding la uncertainties (blue band). The current SAID solution (WI08) [23] 
is represented by the dashed curve; their five single-energy values (at T = 20, 30, 47, 
66, and 90 MeV) are also shown. To enable the meaningful comparison of the values 
contained in this figure, an energy-dependent baseline Sr (= (0.20-T + 1.54)T-10~^, 
with T in MeV and Sfi{T) in degrees) was subtracted from all data. 
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Fig. 4. The em-modified hadronic phase shift 5ii_ (P31) from the present work, 
along with la uncertainties (yellow band). Shown is also the ZUAS12 solution [1], 
obtained on the basis of the elastic-scattering data below 100 MeV, along with the 
corresponding la uncertainties (blue band). The current SAID solution (WI08) [23] 
is represented by the dashed curve. 
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Fig. 5. The em-modified hadronic phase shift d^^ (P13) from the present work, 

along with la uncertainties (yellow band). Shown is also the ZUAS12 solution [1], 

obtained on the basis of the elastic-scattering data below 100 MeV, along with the 

corresponding la uncertainties (blue band). The current SAID solution (WI08) [23] 

is represented by the dashed curve. 
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Fig. 6. The em-modified hadronic phase shift 6i_ (Pll) from the present work, 
along with la uncertainties (yellow band). Shown is also the ZUAS12 solution [1], 
obtained on the basis of the elastic-scattering data below 100 MeV, along with the 
corresponding la uncertainties (blue band). The current SAID solution (WI08) [23] 
is represented by the dashed curve; their five single-energy values (at T = 20, 30, 
47, 66, and 90 MeV) are also shown. 
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Fig. 7. The scale factors Zj for the vr^p data, obtained from the common fit to 
the truncated vr+p and CX databases using the ETH model (see Section 4.2). The 
values, corresponding to the two data sets which were freely floated (see Table 1 of 
Ref. [1]), have not been included. 
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Fig. 8. The scale factors Zj for the CX data, obtamed from the common fit to 
the truncated vr+p and CX databases using the ETH model (see Section 4.2). The 
values, corresponding to the four data sets which were freely floated (see Table 2), 
have not been included. The data sets with the largest scale factors Zj are the three 
FITZGERALD86 data sets which were not freely floated, as well as the FRLEZ98 
data set. 
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Fig. 9. The distribution of the normahsed residuals, obtained from the common fit 
to the truncated tt~^p and CX databases using the ETH model (see Section 4.2). 
Also shown (solid curve) is the optimal Gaussian fit to the data. 
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Fig. 10. The scale factors Zj for free floating (evaluated with Eq.(5) of Ref. [1]) for the 
CX data sets, obtained on the basis of the ZUAS12 solution [1], plotted separately 
for differential cross sections (DCS), total cross sections (TCS), analysing powers 
(AP), and the results for the coefficients of the Legendre expansion of the DCS 
(LEC). The four FITZGERALD86 data sets, which had been freely floated (see 
Table 2), are not shown. Not shown in the figure is also the result of Ref. [17] for 
the isovector scattering length 6i. 
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Fig. 11. Our two predictions corresponding to the CX DCS for centre-of-mass scat- 
tering angle 9 = 0° around the s- and p-wave interference minimum. The ZUAS12 
prediction [1] is represented by the blue band, whereas the ZUAS12a one (this work) 
by the yellow band. Both bands indicate la uncertainties. 
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Fig. 12. Our two predictions corresponding to the CX DCS at 30 MeV; 9 denotes 
the centre-of-mass scattering angle. The ZUAS12 prediction [1] is represented by 
the blue band, whereas the ZUAS12a one (this work) by the yellow band. Both 
bands indicate la uncertainties. 
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Fig. 13. Our two predictions corresponding to the CX DCS at 80 MeV; 9 denotes 
the centre-of-mass scattering angle. The ZUAS12 prediction [1] is represented by 
the blue band, whereas the ZUAS12a one (this work) by the yellow band. Both 
bands indicate la uncertainties. 
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Fig. 14. Our two predictions corresponding to the CX DCS in the kineniatical 
region of the JIA08 [20] experiment. The ZUAS12 prediction [1] is represented by 
the blue band, whereas the ZUAS12a one (this work) by the yellow band. Both 
bands indicate la uncertainties. Not included in the uncertainty of their data is the 
10% normalisation uncertainty of the experiment of Ref. [20]. 



41 



IE 



(a) 



n>Ao 



TC'^ 



n 




% 





Fig. 15. Influence of the ij — n^ mixing on the Feynman diagrams of the ETH model 
in case of the CX reaction. 
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